Mon. Not. R. Astron. Soc. 000, [Tl{9] (2007) Printed 1 February 2008 (MN WT^K. style file v2.2) 



A double main sequence turn-off in the rich star cluster 
NGC 1846 in the Large Magellanic Cloud 
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ABSTRACT 

We report on HST/ACS photometry of the rich intermediate-age star cluster NGC 
1846 in the Large Magellanic Cloud, which clearly reveals the presence of a dou- 
ble main sequence turn-off in this object. Despite this, the main sequence, sub-giant 
branch, and red giant branch are all narrow and well-defined, and the red clump is 
compact. We examine the spatial distribution of turn-off stars and demonstrate that 
all belong to NGC 1846 rather than to any field star population. In addition, the 
spatial distributions of the two sets of turn-off stars may exhibit different central con- 
centrations and some asymmetries. By fitting isochrones, we show that the properties 
of the colour-magnitude diagram can be explained if there are two stellar populations 
of equivalent metal abundance in NGC 1846, differing in age by w 300 Myr. The ab- 
solute ages of the two populations are ~ 1.9 and ~ 2.2 Gyr, although there may be 
a systematic error of up to ±0.4 Gyr in these values. The metal abundance inferred 
from isochrone fitting is [M/H] « —0.40, consistent with spectroscopic measurements 
of [Fe/H]. We propose that the observed properties of NGC 1846 can be explained if 
this object originated via the tidal capture of two star clusters formed separately in a 
star cluster group in a single giant molecular cloud. This scenario accounts naturally 
for the age difference and uniform metallicity of the two member populations, as well 
as the differences in their spatial distributions. 
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1 INTRODUCTION 

The Large Magellanic Cloud (LMC) possesses an extensive 
system of massive star clusters, covering the full range of 
ages 10 6 — 10 10 yr. This age spread, which is not present in 
the Galactic globular cluster population, combined with the 
close proximity of the LMC system, means that these ob- 
jects have proved vital to our understanding of star cluster 
formation and evolution. Furthermore, the age and metal- 
licity distributions of the LMC clusters, together with their 
kinematics, offer important insights into the formation and 
subsequent development of the LMC itself. 

As part of a study of the structural evolution of mas- 
sive stellar clusters, we have conducted a snapshot imaging 
survey of some ~ 50 such objects in the LMC and SMC us- 
ing the Advanced Camera for Surveys (ACS) on-board the 
Hubble Space Telescope (HST). Many of the target clusters 
have not previously been investigated in any detail, and we 
are therefore working on deriving accurate photometric ages 
and metallicities for the sample IjMackev fc Gilmord 120041 : 
iMackev. Payne fc Gilm"orell2006l ). While reducing our obser- 
vations of the poorly-studied rich intermediate-age LMC 
cluster NGC 1846, we recently noticed that the colour- 
magnitude diagram (CMD) for this object exhibited a very 



peculiar main sequence turn-off region; further scrutiny re- 
vealed clearly the presence of two distinct turn-offs. With 
a few notable exceptions (for example, the massive Galac- 
tic globular cluster ui Centauri), the vast majority of star 
clusters which have been resolved by observations into in- 
dividual stars are comprised of single stellar populations - 
that is, groups of stars formed at the same time and with 
the same chemical composition. NGC 1846 is therefore an 
important object, as it can potentially offer us new insights 
into the processes of star and star cluster formation. 

In this paper we present our photometry of NGC 1846 
together with a detailed investigation of its extremely un- 
usual CMD (Sections [2] and I3.1|l . We examine the spatial 
distribution of stars in the main sequence turn-off region, 
and demonstrate that both turn-offs are clearly associated 
with the cluster rather than being an artifact introduced by 
field star contamination (Section 13. 2\ . In addition, we show 
that neither differential reddening, nor a significant line-of- 
sight depth to NGC 1846 can explain the double turn-off. 
From isochrone fitting, we determine the ages and metallic- 
ities of the two populations present in NGC 1846 (Section 
I3.3[l . and use these results to constrain the possible means 
by which this object can have been formed (Section [4}. 
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Figure 1. Drizzled ACS/WFC F814W image of the core of NGC 1846 (exposure 200s). The full image has been cropped to include only 
WFC chip 1 and part of WFC chip 2 — these cover the central cluster region. Note that the brightest stars in the cluster are saturated. 



2 OBSERVATIONS AND DATA REDUCTION 

Our observations were made during HST Cycle 12 on 2003 
October 08 using the ACS Wide Field Channel (WFC). 
As a snapshot target, one frame was taken in each of two 
filters - F555W (dataset j8ne55z9q) and F814W (dataset 
j8ne55zeq). Exposure durations were 300s and 200s, respec- 
tively. The ACS WFC consists of two 2048 x 4096 pixel 
CCDs separated by a gap ~ 50 pixels wide. The plate scale 
is 0.05 arcsec per pixel, resulting in a total areal coverage 
of approximately 202 x 202 arcsec. The core of NGC 1846 
was positioned at the centre of chip 1 so that the inter-chip 
gap did not impact on the innermost region of the cluster. 
In order to help with the identification and removal of hot- 
pixels and cosmic rays, the F814W image was offset from the 
F555W image by ~ 2 pixels in both the x- and y-directions. 

The data products produced by the STScI reduction 
pipeline, which we retrieved via the public archive, have had 
bias and dark-current frames subtracted and are divided by 
a flat-field image. In addition, known hot-pixels and other 
defects are masked, and the photometric keywords in the 
image headers are calculated. We also obtained distortion- 
corrected (drizzled) images from the archive, produced using 
the pyraf task multidrizzle. Part of the drizzled F814W 
image is displayed in Fig. [1] 

W e used the dolphot photometry software (|Dolphinl 
2000), specifically the ACS module, to photometer our 
flatfielded F555W and F814W images, dolphot performs 
point-spread function fitting using PSFs especially tailored 



to the ACS camera. Before performing the photometry, we 
first prepared the images using the dolphot packages AC- 
SMASK and SPLITGROUPS. Respectively, these two packages 
apply the image defect mask and then split the multi-image 
STScI FITS files into a single FITS file per chip. We then 
used the main dolphot routine to simultaneously make 
photometric measurements on the pre-processed images, rel- 
ative to the coordinate system of the drizzled F814W image. 
We chose to fit the sky locally around each detected source 
(important due to the crowded nature of the target), and 
keep only objects with a signal greater than 10 times the 
standard deviation of the background. The output photom- 
et ry from dolphot is on the calibrated VEGAMAG scale 
of H lrianm et all (|2005h . and corrected for charge-transfer 
efficiency (CTE) degradation. 

To obtain a clean list of stellar detections with high 
quality photometry, we applied a filter employing the sharp- 
ness and "crowding" parameters calculated by dolphot. 
The sharpness is a measure of the broadness of a detected 
object relative to the PSF - for a perfectly-fit star this pa- 
rameter is zero, while it is negative for an object which is 
too sharp (perhaps a cosmic-ray) and positive for an object 
which is too broad (say, a background galaxy). The crowd- 
ing parameter measures how much brighter a detected ob- 
ject would have been measured had nearby objects not been 
fit simultaneously. We selected only objects with —0.15 ^ 
sharpness ^ 0.15 in both frames, and crowding ^ 0.5 
mag in both frames. We also only kept objects classified 
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Figure 2. Colour-magnitude diagram for NGC 1846. All detected 
sources which passed successfully through our quality filter are 
plotted (34534 objects). The main cluster sequences are clearly 
visible; however there is some field star contamination present in 
key regions, particularly along the RGB and main sequence. Even 
so, the double main sequence turn-off is strikingly evident. 

by dolphot as good stars (object type 1), as opposed to 
elongated or extended objects (object types > 1). 



3 RESULTS 

3.1 Colour-magnitude diagram 

The colour-magnitude diagram (CMD) for NGC 1846 is pre- 
sented in Fig. [5] All detected stars which passed successfully 
through the quality filter are plotted. There is some evident 
field-star contamination present in this CMD, so in Fig.|3]we 
plot only stars within 30" of the cluster centre, which lies 
near pixel coordinates (x c ,y c ) — (2120,3150) (see Section 
13.21 below). Since this selected region is very much smaller 
in area than the full ACS coverage, field-star contamination 
in the CMD is greatly reduced compared to Fig. [5] and al- 
lows one to clearly identify the primary cluster sequences. 
In principle we could apply a more sophisticated statisti- 
cal subtraction of field stars; however this is unnecessary to 
achieve the goals of the present work. Selecting the inner- 
most cluster regions, as in Fig. [3] clearly offers a sufficiently 
uncontaminated CMD. 

The most striking and unusual feature of the CMD for 
NGC 1846 is the presence of an apparently double main 
sequence turn-off. This is clearly visible in both Figs. [2] 
and O suggesting it is not an artifact of field star contam- 
ination (we address this issue more completely in Section 
13.21 below) . Apart from this, the CMD is as expected for 
an intermediate-age Magellanic Cloud cluster. There is a 
relatively narrow main sequence, a down-sloping sub-giant 
branch (SGB), and a narrow red giant branch (RGB), as 
well as a compact red clump (RC) near mpsssw ~ 19-1 and 



Figure 3. Colour-magnitude diagram for NGC 1846, now plotted 
using only quality- filtered stars at projected radial distances less 
than 30" from the centre (6599 objects). The majority of field star 
contamination has been removed leaving the cluster sequences 
clearly visible, particularly the double main sequence turn-off. 
Note that the view in this plot is focussed on the turn-off region, 
SGB and RGB, so the scale is different from that of Fig. [2] 

"1F555W — wiF8i4W ~ 1-0. On the RGB just brighter than 
the RC level a clear knot of stars is visible, especially in 
Figure [2] This is possibly the RGB Bump - the point at 
which the hydrogen-burning shell in a red giant star reaches 
the discontinuity left in the hydrogen abundance profile due 
to inner penetration of convection. A similar feature has 
recently been observe d in the rich intermedi ate-age LMC 
cluster NGC 1978 by iMucciarelli et al.l l|2007l ) - an object 
which presents a very similar CMD to that for NGC 1846 
(although without the double main sequence turn-off). 

It is worth noting that a third turn-off is clearly visi- 
ble in Fig. [2] near ttifsssw ~ 22.5. This is clearly part of 
the contaminating field star population (since it disappears 
in Fig. [3}, and is representative of a much older population 
than that present in NGC 1846. IMucciarelli et al.1 (|2007[ ) 
observe a similar population in the field surrounding NGC 
1978, and interpret it as the signature of a major star for- 
mation episode which occurred ~ 5 — 6 Gyr ago when the 
LMC and SMC were gravitationally bound. 

Above mp 5 55w ~ 18.0, the width of the cluster RGB 
appears to increase significantly. Examining the data qual- 
ity flags produced by dolphot we found that all stars above 
this level were flagged as being saturated in their central pix- 
els (see e.g., the brightest stars in Fig.[T|. The accuracy of all 
photometry on the RGB above rriFsssw = 18.0 is therefore 
dubious, and the apparent increased spread is most likely 
artificial. 

Based on an unpublished CMD. lGrocholski et all (|2006T ) 
suggested that NGC 1846 suffers from differential reddening; 
however, from the narrow (lower) RGB and the compact RC 
visible in Fig. [3] we see no evidence for significant differential 
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Figure 4. Colour-magnitude diagram of the main sequence turn- 
off region of NGC 1846, with all stars plotted as in Fig. [2] The 
two marked boxes were determined empirically using the CMD 
in Fig. \3\ and are designed to isolate the "upper" and "lower" 
turn-off stars. The main sequence itself has been excluded to try 
and maintain as clean a separation as possible. 



reddening, at least in the central regions of the cluster. Cer- 
tainly the apparent double main sequence turn-off cannot 
be produced by such an effect without introducing a large 
spread into the other features on the cluster CMD. Simi- 
larly, the observed lack of such a spread in these features 
militates against a large line-of-sight depth in the system 
being responsible for the properties of the turn-off region. 



3.2 Spatial distribution of turn-off stars 

The most intriguing aspect of the CMD for NGC 1846 is 
the double main sequence turn-off. Such a clearly defined 
example has not previously been observed in a Magellanic 
Cloud cluster - indeed, it is rare to find such a feature in any 
type of star cluster. The turn-off region is suggestive of two 
separate cluster populations. If this is so then it is important 
to ascertain whether there are any differences in their spatial 
distributions, as is observed, for example, for the multiple 
populations found in the peculiar Galactic globular cluster 
uj Centauri (see e.g.. IVillanova et alj|2007l . and references 
therein). 

We first need to demonstrate that all stars in the main 
sequence turn-off region belong to the cluster. To this end, 
we used Fig. [3] to isolate the areas on the CMD covered by 
the two visible turn-offs. We then searched through the full 
list of stars plotted in Fig. [2] and chose only those lying in 
these two regions. The results of this process are displayed 
in Fig. U The two boxes define the sets of "upper" and 
"lower" turn-off stars. By examining the spatial distribu- 
tions of these two ensembles, we can check that they both 
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Figure 5. Upper panel: Pixel coordinates of all 34534 detected 
stars with high quality photometry, along with our measured cen- 
tre at (x c ,y c ) = (2120, 3150) (solid point). The solid circle marks 
a radius of 80" (1600 pixels) about this centre. Lower panel: Dis- 
tributions of projected radial distance from the cluster centre for 
all stars (solid line), upper turn-off stars (dashed line), and lower 
turn-off stars (dot-dashed line). The vertical dotted line marks a 
radius of 80", as plotted in the upper panel. Approximately 90 
per cent of the turn-off stars lie within this radius. 



do belong to the cluster, and search for differences between 
them. 

In Fig. [5] we plot the distribution of projected radial 
distance from the cluster centre for all stars and com- 
pare this with the distributions of projected radial distance 
from the cluster centre for the upper and lower turn-off 
stars. We determined the cluster centre using a surface- 
bri ghtness testing algorithm very similar to that described 
bv lMackev fc Gilmorei (|2003h , with the result that the clus- 
ter centre lies within ~ 20 pixels (1") of the pixel coordi- 
nates (x c ,y c ) = (2120,3150). The pixel coordinates of all 
stars are plotted in the upper panel of Fig. [5] along with a 
point marking the measured cluster centre. 

The radial distribution of all stars (including both clus- 
ter and field stars) , plotted as a solid line in the lower panel 
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Figure 6. Contour plots of the spatial distribution of upper turn-off stars (top left panel), lower turn-off stars (top right panel), and 
all detected stars (lower panel). Axis units are pixel coordinates - the plots are zoomed on the cluster centre, so do not cover the entire 
ACS field of view. They were produced by calculating a 2D histogram of stellar positions using a box size of 300 X 300 pixels, and then 
smoothing and contouring the result. The marked contours in the upper two panels delineate approximately equivalent density levels. 



of Fig. [5] is much less centrally concentrated than those of 
the upper and lower turn-off stars (dashed line, and dot- 
dashed line, respectively). Approximately 90 per cent of the 
turn-off stars lie within a radius of ~ 80" (1600 pixels) of the 
cluster centre, as marked in the upper panel of the Figure. 
This observation strongly implies that these objects are as- 
sociated with the cluster itself rather than representing any 



type of field star contamination. There is also apparently a 
small difference between the radial distributions of the two 
sets of turn-off stars, suggesting that the lower turn-off stars 
may be a little less centrally concentrated than the upper 
turn-off stars. 

In Fig. [6] we show contour plots of the spatial distribu- 
tions of all stars (lower panel) and of both sets of turn-off 
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stars (upper two panels). This allows us to further compare 
the spatial distributions of the turn-off stars to each other 
and to the rest of the cluster stars. We produced each plot 
by first making a 2D histogram of the pixel coordinates of 
the stars in the relevant group, using a bin size of 300 x 300 
pixels. Selection of this size represents a trade-off between 
high spatial resolution and having a significant number of 
stars per bin for the turn-off plots; we retained the same 
box size for the plot including all stars to allow direct com- 
parison. After calculating the histogram we smoothed the 
result using the MIN_CURVE_SURF function in idl, and then 
obtained a contour map using the CONTOUR procedure in 
idl. The results still inevitably contain a small amount of 
residual boxiness from the 2D histogram; however this does 
not affect the overall shapes and spacing of the displayed 
contours. 

The lower panel of Fig. [6] shows that the cluster con- 
sidered as a whole appears to be relatively smoothly and 
symmetrically distributed. The centre matches well our mea- 
sured value of (x c , y c ) = (2120, 3150), and there is no strong 
evidence for significant ellipticity or dumpiness. From the 
upper two panels, it is clear that the two sets of turn-off 
stars are very definitely associated with the cluster, as we 
surmised from Fig. [S] In addition, within the accuracy of 
centroiding from the maps (roughly ±50 pixels, or ±2.5") 
their centres are coincident with each other and the cluster 
as a whole. It is possible that the set of upper turn-off stars is 
more centrally concentrated than both the set of lower turn- 
off stars and the cluster as a whole (this would be consistent 
with Fig.[SJ), and furthermore, that the distribution of upper 
turn-off stars may also exhibit some significant asymmetry. 
However, both improved photometry of the turn-off region 
and a more sophisticated statistical analysis are required to 
draw a firm conclusion on these issues. 

Together Figs. [5] and [6] offer strong evidence that both 
sets of turn-off stars are definitely associated with NGC 1846 
rather than any field star population. Furthermore, there 
are tantalizing indications that the spatial distributions of 
the two sets of turn-off stars may be somewhat different to 
each other, although improved measurements are required 
to confirm this. 



3.3 Isochrone fitting 

From the above analysis, we conclude that the observed 
double main sequence turn-off in NGC 1846 results from 
the presence of two separate stellar populations in this clus- 
ter. However, even though the two main sequence turn-offs 
are quite distinct, the RGB (below the saturation level at 
m-F555w = 18.0) is narrow, and the RC is of a compact na- 
ture. This implies that the two populations cannot be of 
strongly differing metallicities, since both the RGB and the 
RC of a more metal rich population would lie considerably 
further to the red than those for a more metal poor popu- 
lation, and the RC would also be somewhat less luminous 
for the more metal rich population than for the metal poor 
population. 

To quantify the main characteristics (age, metallicity) 
of the two populations in NGC 1846, we fit isochrones 
to the cluster CMD. We used two sets of stellar mod- 
els which have isochrones calculated in the F555W and 
F814W ACS/WFC filter systems - those of the Padova 



group JCirardi et alfe oOO). and those from the BaSTI evolu- 



tionary ^ode^^^^^^^MoO^^e^^^MQM)- We 

selected the "basic set" of Padova tracks with solar-scaled 
distribution of metals (in particular, that is, not enhanced 
in a-elements), which inclu de some degree of convective 
overshooting (for details see iGirardi et all l2000h . and the 
"non-canonical" solar-scaled BaS TI tracks, which also e m- 
ploy convective overshooting (see IPietrinferni et al.ll2004l ). 

Apart from the double turn-off, NGC 1846 presents a 
rather similar CMD to that of NGC 19 78, which was re- 
centl y measured to be 1.9 ± 0.1 Gyr old (Mucciarcll i et all 
I2007T ). There are two available spectroscopi c measurements 
of the metallicity of NGC 1846 - that of lOlszewski et al.1 
|l99ll ) who found [Fe/H] = -0.70 ± 0.20 from one RGB 
star, and that of iGrocholski et all |2006) who obtained 
[Fe/H] = -0.49 ± 0.03 from 17 RGB stars. For the Padova 
models we used the interactive web form available on the 
group's web site to construct a fine grid of isochrones about 
these values, sampling an age range 1.0 ^ r 3.0 Gyr at in- 
tervals of 0.1 Gyr, and a metal abundance range 0.00250 ^ 
Z sC 0.00950 at intervals of 0.00025. The total metallicity 
[M/H] = log(Z/Z©) where Z Q « 0.019, so this abundance 
range corresponds to —0.88 ^ [M/H] ^ —0.30. Assuming 
the a-element enhancement in NGC 1846 is small, which 
would be con sistent with other inter mediate-age LMC clus- 
ters (see e.g. [MucciarenrelaD ioo3, for NGC 1978), then 
[M/H] ~ [Fe/H]. For the BaSTI tracks, we used the web 
form available on that group's web site to construct a grid 
of isochrones sampling the same age range as the Padova 
grid, at the same intervals, but for only two metal abun- 
dances - Z = 0.004 and Z = 0.008 (the BaSTI web form 
does not allow interpolation in metal abundance). 

We first attempted to find the best fitting isochrone for 
the upper main sequence turn-off. We did this by locating 
by eye three fiducial points on the CMD of the cluster's 
central region: the magnitude and colour of the turn-off, 
the magnitude of the RC, and the colour of the RGB at 
a level 1.0 mag brighter than the level of the turn-off. This 
latter point was selected simply as a point lying on the lower 
RGB at a level intermediate between that of the red end 
of the SGB and that of the RC. We then calculated the 
intervals Av and Ac, which are, respectively, the difference 
in magnitude between the level of the turn-off and the level 
of the RC, and the difference in colour between the turn- 
off and the RGB fiducial point. These intervals are useful 
because Av is strongly sensitive to cluster age (and weakly 
to cluster metallicity), while Ac is sensitive to both cluster 
age and metallicity. We obtained Av = 1-65 ± 0.1 mag and 
A c = 0.68 ± 0.01 mag. 

Next, we calculated the same intervals for all isochrones 
on each of our two grids, and selected only those with values 
lying within certain tolerances of the cluster measurements 
(±0.2 mag for Av and ±0.02 mag for Ac). This resulted in 
lists of fewer than ten isochrones per grid, which we then fit 
to the CMD by eye. To do this, we calculated the offsets in 
magnitude and colour required to align the turn-off of the 
isochrone with that of the CMD, and then those required 
to align the RC of the isochrone with that of the CMD. 
This resulted in two offsets in magnitude, and two offsets 
in colour, which we averaged and then applied to overplot 
the isochrone on the CMD. Given the narrowness of most of 
the sequences on the CMD, it was then straightforward to 
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Figure 7. Best-fitting Padova isochrones overplotted on the clus- 
ter CMD. The upper and lower turn-off isochrones have the same 
metal abundance Z = 0.0075, and ages t u = 1.5 Gyr and t; = 1.8 
Gyr, respectively. Both have been shifted by E(V — I) = 0.045 
and /J, = 18.37 to achieve the best-fitting alignment. 



identify the best fitting isochrone for each grid. The resulting 
offsets Sy and Sc provide estimates for the distance modulus 
to the cluster (n) and the foreground extinction, in that 
5c = J?(mF555w — ntF 8i4w) ~ E(V — I) (sin ce the extinction 
values are small, see ISirianni et alj 120051 ) and Sv ~ fJ. ± 
2.37E(V-I). 

We then moved to the lower turn-off and repeated the 
above process. In Section 13.21 we demonstrated that both 
populations in NGC 1846 share an approximately common 
spatial centre. Furthermore, the narrowness of the CMD se- 
quences implies that there is little or no differential redden- 
ing towards the cluster centre, nor a significant line-of-sight 
depth to the system. In aligning the best-fitting isochrones 
to the lower turn-off, we therefore applied the same distance 
modulus and extinction values calculated for the upper turn- 
off. 

The best-fitting Padova isochrones may be seen in 
Fig. [7] We found that the highest quality solution was for 
isochrones of the same metal abundance (Z = 0.0075, or 
[M/H] = -0.40) but differing ages: r u = 1.5 Gyr for the 
upper turn-off and t; = 1.8 Gyr for the lower turn-off. The 
isochrone filtering process outlined above provided an in- 
dication of the random uncertainties associated with these 
results - isochrones with ±0.00025 in Z and ±0.1 Gyr in 
age gave reasonable fits, but clearly not as good as for the 
selected values. The fact that the two isochrones have the 
same metallicity is primarily driven by the narrowness of the 
observed RGB and the compact nature of the RC. Matching 
these features is successfully achieved if there is solely an age 
difference between the two populations. The observed width 
of the SGB is successfully matched at the red end, but is 
apparently a little too wide at the blue end. The reason for 



this is not clear - it may represent a shortcoming of the stel- 
lar models (although we note that the Basti isochrones also 
predict too wide an SGB at the blue end), or may represent 
another very unusual feature of the CMD for NGC 1846. 

The distance modulus and colour excess obtained for 
the fit described above are fi = 18.37 and E(V - 1) = 0.045. 
These are both slightly lower than the ca nonical values for 
the L MC; a similar effect was noted by iMucciarelli et all 
(2007) when fitting Padova isochrones to their CMD for 
NGC 1978. 

Fig.|8]shows the best-fitting BaSTI isochrones. For these 
models we did not have the luxury of fine sampling in metal 
abundance; however, Z — 0.008 is close to the best-fitting 
Padova metallicity, and provides adequate results. The best- 
fitting BaSTI isochrones are somewhat older than those for 
the Padova models - we found r u = 2.2 Gyr for the up- 
per turn-off and n — 2.5 Gyr for the lower turn-off. Un- 
certainties are similar to those for the Padova model fit- 
ting. The RGB, RC, and width of the SGB are slightly 
more successfully matched than with the best-fitting Padova 
isochrones; however, the turn-off regions are not quite as well 
matched. This may be an artifact of having a slightly incor- 
rect metallicity for the BaSTI models. The distance mod- 
ulus and colour excess obtained for the fit are /j, = 18.57 
and E(V — I) = 0.075. In this case, the distance modulus is 
slightly higher than the canonical LMC value. 

While our two pairs of measured absolute ages show 
an offset of ~ 0.7 Gyr, both sets of models show that 
the appearance of the NGC 1846 CMD can successfully 
be reproduced by adopting two stellar populations of the 
same metallicity and a difference in age of ~ 300 Myr. 
Our derived metallicity [M/H] ~ —0.40 is consistent with 
the spectroscopi c iron abundance measured recently by 
iGrocholski et all (|2006l ) from a sample of 17 RGB stars. 



4 DISCUSSION AND CONCLUSIONS 

Our HST/ACS photometry of the rich intermediate-age 
LMC cluster NGC 1846 has revealed an unusual CMD, ex- 
hibiting two main sequence turn-offs, but otherwise narrow 
and clearly defined primary features (RGB, SGB, RC). By 
examining the spatial distributions of the two sets of turn-off 
stars, we have shown that both are unequivocally associated 
with the cluster rather than any field-star population. The 
set of upper turn-off stars is possibly more centrally concen- 
trated than that of the set of lower turn-off stars, and may 
be asymmetrically distributed. Even so, the two sets of turn- 
off stars share an approximately common centre. Via careful 
fitting of isochrones from two different groups of stellar evo- 
lution models (Padova and BaSTI), we have demonstrated 
that the observed features on the CMD can be explained 
by assuming that NGC 1846 harbours two stellar popula- 
tions of the same metal abundance ([M/H] ~ —0.40) but 
differing in age by ~ 300 Myr. The absolute ages of the two 
populations are t u — 1.9 ± 0.1 Gyr and r; = 2.2 ± 0.1 Gyr. 
These are a straight average of the ages from the two differ- 
ent isochrone sets, and the uncertainties represent the ap- 
proximate random uncertainties associated with the fitting 
process. We note that additional, larger systematic uncer- 
tainties are likely present - the two pairs of absolute ages 
are offset by 0.7 Gyr between the two sets of isochrones. 
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Figure 8. Best-fitting BaSTI isochrones overplotted on the clus- 
ter CMD. The upper and lower turn-off isochrones have the same 
metal abundance Z = 0.008, and ages t u = 2.2 Gyr and t; = 2.5 
Gyr, respectively. Both have been shifted by E(V — I) = 0.075 
and fx = 18.57 to achieve the best-fitting alignment. 



Given these results, it is natural to ask how such an 
object as NGC 1846 might have been formed. We can think 
of two options - either this cluster underwent two distinct 
episodes of star formation separated by 300 Myr, or it is 
the result of the merger of two star clusters formed 300 Myr 
apart. 

Regarding the possibility that NGC 1846 underwent 
two separate episodes of star formation, it is difficult to 
see how it can have retained enough gas for this to hap- 
pen. There is strong evidence that the combined effects of 
massive stellar winds and supernova explosions expel any 
remaining gas from a very young star cl uster within only a 
few Myr (e.g.. iBastian fe Goodwm]|2006l ). leading to strong 
chan ges in the cluster's internal dynamics and structure (see 
e.g., iGoodwin fe Bastianl 120061 . and references therein). A 
cluster must be very massive to retain any gas in its po- 
tential well. The only clusters where multiple episodes of 
star formation appear evident are globular clusters at the 
very upper end of t he globular cluster m ass function, such 
as u Centauri (e.g., rVillanova et al.ll2007h . and Gl in M31 
jMevlan et al.ll200ll ). NGC 1846 does not fulfil this crite- 
rion. Furthermore, in all such cases where multiple episodes 
of star formation have taken place, self-enrichment has ap- 
parently occurred. The timescale for such enrichment is not 
well defined; however, it seems possible that 300 Myr could 
be long enough given that significant chemical processing 
must have already occurred in the most massive cluster stars 
(in particular in type II supernovae and massive asymptotic 
giant branch stars), within the first ~ 100 Myr. It is there- 
fore not clear why the two populations in NGC 1846 should 
apparently have the same metal abundance. 

The merger scenario fares better. It is well known 



that both the LMC and SMC possess populations 
of candidate binary (a n d mu lti ple) star clusters (e.g., 
Bhatia & Hatzidimitrio u| Il988l ; lHatzidimitriou fc Bhatia! 
Il990h . and by using statistical arguments it can be shown 
that a signific ant fraction of these are likel y to be phys- 
ically linked l|Bhatia fc Hatzidimitrioul Il988l ). Gravitation- 
ally bound star clusters are predicted to merge on time- 
scales of tens of Myr to a few hundred Myr (e.g., 
iPortegies Zwart fc Rushll2006l ). which is roughly consistent 
with the observed age difference between the two popula- 
tions in NGC 1846. Old binary pairs - such as SL 349-SL 
353 with an age ~ 500 Myr - are certai nly observed in 
the LMC l|Leon. Bergond fc Valrenarilll999h . However, the- 
ories for the formation o f bound binary star clusters (e.g., 
iFuiimoto fc Kumail 11997m suggest that such objects should 
be roughly coeval, with age differences less than ~ 100 Myr 
at most. This is inconsistent with the observed populations 
in NGC 1846, and indeed with several known binary pairs 
in the LMC - such as SL 356-SL 357, which have ages of 70 
Myr and 600 Myr, respectively (|Leon et al.lll999l '). 

One way to circumvent this is if some clusters are 
formed in star cluster groups (SCGs) in giant molecular 
clouds (see e.g. jLeon et al.ll"l999l . and references therein). In 
such complexes, cluster formation can be spread over several 
hundred Myr. Furthermore, clusters observed to be physi- 
cally associated at present do not have to have been formed 
this way, as the cross-section for tidal capture is much in- 
creased due to the relatively high density of star clusters 
within the SCG. This scenario therefore allows for the exis- 
tence of cluster pairs which are older than the typical time- 
scale for merging, as well as for cluster pairs with large age 
differences, both of which types are seen in the LMC. If the 
two components of NGC 1846 were formed in a SCG in the 
same giant molecular cloud, this also naturally explains their 
matching metal abundances. 

In this scenario, the two clusters comprising the present 
day NGC 1846 must have merged less than 1.5 — 2.2 Gyr 
ago, based on our derived ages for the upper turn-off popula- 
tion. The median relaxation time in intermed iate-age LMC 
clusters is typically of order 1 — 2 Gyr (e.g., iMackev et al.l 
120071 ); the central relaxation time may be up to a factor 
ten shorter than this. The very inner region of NGC 1846 
is therefore quite possibly dynamically old enough for two 
merging clusters to now be sharing a common centre and 
be fairly well mixed through relaxation processes; however 
it is probable that asymmetries could still exist outside the 
core, as our observations have hinted. It would be extremely 
interesting to investigate the evolution of a merged LMC- 
type cluster via realistic TV-body modelling to examine if 
the observed features of NGC 1846 can be reproduced. Such 
calculations may allow the time when merging occurred to 
be constrained - this would be extremely useful in the con- 
text of examining the SCG scenario. Direct TV-body mod- 
elling of LMC-type clusters is now possible - IMackev et al.l 
|2003) present models with TV = 10 5 particles integrated 
over a Hubble time of evolution. Since we do observe asym- 
metries in the distribution of the two populations in NGC 
1846, dynamical signatures offering clues to the formation 
process may also still be present outside the cluster centre. 
It would certainly be interesting (although observationally 
challenging) to investigate the internal dynamics of this ob- 
ject. TV-body modelling would be helpful in this regard as 
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it would allow the expected level of the signature to be pre- 
dicted. NGC 1846 patently warrants further study, including 
improved photometry of the turn-off region, for important 
insights into star cluster formation processes. 
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